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We calculate the low-frequency magnetoplasmon excitation spectrum for a square array of quan- 
tum dots on a two-dimensional (2D) graphene layer. The confining potential is linear in the distance 
from the center of the quantum dot. The electron eigenstates in a magnetic field and confining po- 
tential are mapped onto a 2D plane of electron-hole pairs in an effective magnetic field without any 
confinement. The tight-binding model for the array of quantum dots leads to a wavefunction with 
inter-dot mixing of the quantum numbers associated with an isolated quantum dot. For chosen con- 
finement, magnetic field, wave vector and frequency, we plot the dispersion equation as a function 
of the period d of the lattice. We obtain those values of d which yield collective plasma excitations. 
For the allowed transitions between the valence and conduction bands in our calculations, we obtain 
C/3 plasmons when d < 100A. 
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I. INTRODUCTION 



A two-dimensional (2D) honeycomb lattice of carbon atoms that form the basic planar structure in graphite 
(graphene) has recently been produced 0, 0. Unusual many-body effects in graphene have been attributed to low- 
lying excitations in the vicinity of the Fermi level0, 0, [H, |(|. The influence of an external magnetic field on the 
many-electron properties of graphene results in an unusual quantum Hall effect. The integer quantum Hall effect 
(IQHE) was discovered in graphene in recent experiments [7|, |8|, [9(. The quantum Hall ferromagnetism in graphene 
has been studied theoretically [nj. In addition, the influence of magnetic field on the Bose-Einstein condensation 
and superfluidity of indirect magnetoexcitons in graphene bilayer has been studied [TjJ ■ The spectrum of plasmon 
| excitations in a single graphene layer immersed in a material with effective dielectric constant £(, in the absence of 
t— I . magnetic field (B = 0) was calculated in Refs. [1, [H, [12] ■ (See also [T3, [Tj. Hi l [l 6ll.) The collective plasma excitations 
t-H ' in layered graphene structures in high magnetic field were obtained in Ref. [17| where the instability of these modes 
, was investigated. Recently, several works were reported for quantum dots in graphene [H, [2(| ■ The transport 
characteristics of quantum dot devices etched entirely in graphene have been studied experimentally [2~l| . According 
0^ ", t° Ref- [2l|, these quantum dots at large sizes behave as conventional single-electron transistors. This is just one of the 
areas of research in the fast-growing field of the electronic transport , thermal and optical properties of graphene which 
may have important device applications because of its high mobility (22J. The single-particle states and collective 
modes in semiconductor quantum dots have long been a subject of interest to both theoreticians and experimentalists 
[HI, HH HH ■ The spectrum of collective plasma modes in a 2D array of quantum dots in semiconductors has been 
calculated [26| . 

In this paper, we obtain the dispersion relation for magnetoplasmons in a square array of quantum dots formed 
by a confining potential which is a linear function of distance from the center of the quantum dot. A perpendicular 
magnetic field B is applied. The conditions for the existence of the collective excitations will be investigated. 

This paper is organized as follows. In Sec. [Til we introduce the quasiparticle electron-hole representation for the 
eigenvalue problem for an electron in a single graphene quantum dot in magnetic field. In Sec. IIIII we represent the 
calculations of the magnetoplasmon spectrum in a quantum dot array in graphene. A brief discussion of the results 
of our calculations plasmon instabilities in graphene is presented in Sec. IIV1 
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II. SINGLE-ELECTRON STATES FOR A GRAPHENE QUANTUM DOT IN MAGNETIC FIELD 



We now consider an electron in a graphene layer in the presence of a perpendicular magnetic field B and a confining 
potential U(r') = Ao|r'|, where Ao > is the slope of the potential. The effective- mass Hamiltonian of an electron 
in the absence of scatterers in one valley in graphene located in the x t/-pl ane is given by a 2 x 2 matrix Hamiltonian 
-ff(o) with zero along the diagonal and off-diagonal elements fc x =p in y (27l l28l|. Here, we neglect the Zeeman splitting 
and assume energy degeneracy with respect to the two valleys. Also, in our notation, tt = —iKV + eA, A is the 
vector potential of an electron, vf = V3at/(2h) is the Fermi velocity of electrons with a = 2.566A denoting the 
lattice constant and t 2.71eV is the overlap integral between nearest-neighbor carbon atoms in graphene [29| . 
The eigenvalue problem of an electron in a linear confining potential can be mapped onto one for a noninteracting 
electron-hole pair in an effective magnetic field B c fi under conditions we introduce below. 

The eigenfunction ip T for the Hamiltonian -ff(o) for an electron- hole pair in an effective magnetic field B e g, is also 

the eigenfunction of the magnetic momentum P has the form [30L l3ll . l32| and is given by 



■0 P (R,r) = exp 



e[B cff x r] 



i(r-po) 



(1) 



where R = (r e + r/j)/2, r = r e — r/ t and po = [B c ff x P]/ (ei? 2 ff ). The cylindrical gauge for vector potential is used 
with A e (/j) = f/2[B G ff x r e (7j)]. The wavefunction of the relative coordinate $(r) can be expressed in terms of the 



2D harmonic oscillator eigenfunctions <fr nit n 2 (r). For an electron in Landau level n + and a hole in level n. 
four-component wavefunctions for the relative motion in graphene are [U 



the 



(r) = (V2 



f s+s_<I>| n+ |_ 1> _ hl (r) \_ 

S+*|n+|-l,|n_|(r) 



|n+ , |n_ | — 1 



(r) 



(2) 



\ *|n+|,| n _|(r) / 



where s± = sgn(n±). The corresponding energy of the electron- hole pair E n+ n _ of the Hamiltonian Ht^ is given by 

M 



E,, 



V2 ((hv F /r* B )) 



sgn(n + ) VF+I - sgn(n„ ) 



(3) 



where r* B = \J h/(eB c g) is an effective magnetic length. The 2D harmonic oscillator wavefunctions $m,n 2 ( r ) are given 

by m 



^7, 



,(r) = (2vr)- 1 / 2 2-l m l/ 2 



1 



x — ; — r exp 



V«i!n 2 ! r* B 



sgn(m) r 



2r*, 



(4) 



where is a Laguerre polynomial, m = n% — n,2, n = min(ni, n-i) and sgn(m) m — ► 1 for m = 0. The electron- 

hole eigenfunction given by Eqs. (H|), ©j O along with the energy eigenvalues are the same as the eigenfunction and 
the eigenvalue of the Hamiltonian if(o) for an electron in a quantum dot in external magnetic field. The magnetic 
momentum P = — e[B c ff x r] must be the same for the electron and the electron-hole pair at fixed relative coordinate 
r. The effective magnetic field B c g depends on the slope Ao of the confining potential as well as the external magnetic 
field B and is given by 



B r 



(B 2 + (4A 2 /(e 2 ^)) 



t V2 



(5) 



In the case when there is no electron confinement, i.e., Ao = 0, the effective magnetic field is the same as the applied 
magnetic field and B c s = B. If there is no external magnetic field, i.e., B = 0, then there is still an effective magnetic 
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field due to the electron confinement and is given by B c g = 2Xo/ (evp). The coordinate of the electron-hole relative 
motion r = \r e — 17J is related to the distance r' of an electron from the center of the confining potential by r — r'/2. 

It should be emphasized that the electrostatic potential for a graphene quantum dot was chosen to have a linear 
dependence on the coordinate variable |r'|. This was done because the eigenvalue problem of the Hamiltonian 
described by the Dirac-like equation with a linear potential can be reduced to the Klein- Gordon- type equation with 
a parabolic potential U(r') = Agr' 2 . This Klein-Gordon-type equation can then be mapped onto the electron-hole 
problem in an effective magnetic field _B c ff defined above under the conditions presented. 



III. MAGNETOPLASMONS FOR A QUANTUM DOT ARRAY IN GRAPHENE 

We now turn our attention to an infinite periodic 2D array of quantum dots in a graphene plane, shown schematically 
in Fig. [TJ The quantum dots are formed by a periodic linear confining potential defined by U(r') = ^o| r ' — r j|j 
where Yj is the position vector of a quantum dot. We consider this array with the period d in a perpendicular magnetic 
field B as shown schematically in Fig. [JJ For this system, we apply the tight-binding approximation [26|, [34], [35| (see 
Eq. (2.1) in Ref. [H) 

*k(a?,») = |k,a) =^e lk r ^a(r-r J )exp ^—A-i^J , (6) 
r j 

where k is an in-plane wave vector, ip a is the electron eigenfunction in one quantum dot (which can be represented 
by the wavefunction of an electron-hole pair in an effective magnetic field as ip a — (r/2), where $„ + „_(r/2) 

is given by Eq. @. The index a = {n +1 n_} is a composite quantum number with n + and n_ labeling the electron 
and hole energy levels, respectively. Also, A is the vector potential for the externally applied magnetic field B. 

Note that the electron wave function in periodic systems in a magnetic field is the eigenfunction of both the 
Hamiltonian and the operator of the magnetic translation. The electronic properties in the presence of magnetic field 
are determined by the magnetic flux through a unit cell. [361 ] If this flux (measured in flux quantum (j)Q = h/e) is a 
rational number Q/ifto = p/q = Ba 2 e/2i:h where p and q are prime integers, then the electron wave function which 
is also an eigenfunction of the magnetic translation operator [33] obeys the Bloch-Peierls conditions in the Landau 
gauge which may be expressed as[36j 



* k (a?, y) = * k (x + qa, y + a )e~ lk:cqa e~ ikya e~ 27Tipv / a . (7) 

For integer mi, m%, the vectors &. m (m\qa, n%a) define the magnetic lattice of the crystal. The magnetic Brillouin zone 
is defined by the inequalities —n/qa < k x < ir/qa and 7r/a < k y < k y < n/a. 

Let us demonstrate that the electron wavefunction given by Eq. ^ obeys the Bloch-Peierls conditions in the Landau 
gauge. In this gauge, the vector potential may be chosen as A = (0, Bx, 0). Then, Eq. ([6]) yields 



(8) 



* k (a; + qa,y + a) = ^2e lk ' Tj ip a ((x + qa)i + (y + a)j - r^exp (- l -^-B(x + qa)j ■ . 
Now, let Yj — qai — aj = r'j . Then, Eq. © becomes 

* k (x + qa, y + a) = ^2 e <k ' ( ^ +go?+ °^Va( r ~ r j) ex P (~J~ B ( X + qa ^ ' K + qai + a ^ 
= e lk * qa e lk v a e*' r ^a(^ + VJ-r' j ) 

r' 

3 

(—ie . . „ „ px I r , • j ; \ \ 
——A • r — 2-Kip — 2m 2mp — 
n J q a \ a I I 

= e lfc ^ a e 4fc " Q exp ( -2ir£-) * k (x,y) . (9) 
V 1 a J 



FIG. 1: A periodic 2D array of quantum dots in a graphene plane formed by linear confinement with a spacing d in a 
perpendicular magnetic field B. 

But, this result in Eq. (J9J) is equivalent to Eq. (|7|). The reason is that they coincide when we rotate axes in the frame 
before doing the displacement, i.e., rc-axis — > — y-axis or, replacing x —y/q on the right-hand-side of Eq. ([9]). This 
is now explicitly demonstrated below. Using the gauge A = (—By ,0,0), we obtain 

^(x + qa,y + a) = ^ e^^+^+^V^r - r^) exp (~jB(-y - a)i ■ (r$ + qai + aj)j 

3 

= e ik * qa e lk y a e lk ' r ^«(^ + VJ - r£) 
r 'j 

I —ie . y ~ " 

x exp — — A • r ■ + 2irip + 2irip — h 2ni 



h J ' a ij 1 a 



x e*' r ^a(^ + Vj - r'j) exp |2tu| J | . (10) 

The last line in Eq. (|10[) satisfies the Eq. ||7J) perfectly if we take k to be k modulo 2ir/aq. Therefore, we conclude 
that the electron wavefunction given by Eq. [S] obeys the Bloch-Peierls conditions in the Landau gauge. In Ref. [38[ 
quantum states of a similar system were considered. The energy spectrum in this case consists of magnetic subbands 
defined in a magnetic Brillouin zone. The structure of magnetic subband depends on the number of magnetic flux 
penetrating the unit cell. The wave function given by Eq. ([6]) corresponds to one in Ref. [38| if we do not include 
spin-orbit coupling in our model. 

We do not consider magnetic fields where the flux is p/q in our numerical calculation, since this would involve a 
formalism like the one presented in Refs. [H, H(|. Furthermore, the resulting fractal nature of the energy spectrum 
under these conditions can only be realized numerically for magnetic fields which are several tens of Tesla. This is 
beyond the interest and scope of the present paper. 

Following the procedure adopted in Ref. |26j , we introduce the quantity 



9s9v \ fo(E a ) - fa{E p ) 



where g s — 2 and g v — 2 are the spin and valley degeneracies in graphene, uj denotes frequency, a = (n+,n_), 
(3 = (n',,n'_), E a m\ are the electron-hole pair eigenenergies for the Landau levels in an effective magnetic field. Also, 
fo(E) is the occupation function of the two-particle state. At high temperatures or weak magnetic field, when the 
separation between Landau levels is small so that ksT ^> ?wf/t* b , the occupation of the Landau levels is given by 
the Fermi-Dirac distribution function fo(E a ) = (exp [E a /(kBT)] + 1)~ . 

We now introduce the overlap integral _F a| a(q) involving the eigenstates with labels a = {n+, n_} and (3 — {n' + , n'_} 
through the equation 26] 
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X-F| (0) , i i I -ii/i i i / i n (q) 

+s+s + -F 1 |ij|_i ! | ri _| i |„' + |_i i |„'_|(q) 

+s_s' ^ (0) , . , , Jq) 

' — |n + |,|n_ | — ljn^ | , \n_ \ — 1 v"/ 

+< !|,|n_|,K|,KjW ' ( 12 ) 



where 



ni,n 2 ,nj 



^(q) = | d 2 r e ^ r $;, )n ,(r)$ ni , n2 (r) (13) 
and $„ ljT j 2 (r) is given by Eq. (H|). The angular integral in Eq. (|13p can be carried out using the result 

d9 e ^ cos ( e )-™ e = 2Tr{-i) m J m {-(3) (14) 



with m an integer. If m is even, J m (— (3) = J m (0). 

We employ the procedure introduced in Ref. [261 ] for determining the dispersion relation u> — w(q) for the magneto- 
plasmon excitation frequencies. This can be calculated from the condition of the vanishing of the determinant of the 
matrix with elements CV,-(q, u>). In this notation, the subscript i denotes the pair of quantum numbers a and (3 for 
two different composite electron-hole pairs. Also, j denotes the quantum numbers a' and /?' for two other composite 
electron-hole pairs. The matrix elements dj are defined as 



Cjj(q,w) = Sij 

^e s |q+G| 

In this notation, s s — AneoCb, where eb is the average dielectric constant of the medium where the layer of graphene is 
embedded. We also introduced G = 2ir/d (n x , n y ) which is a reciprocal lattice vector of the square array of quantum 
dots with period d and n x , n y = 0, ±1, ±2, 

In solving the dispersion equation for plasmon excitations, we chose the slope of the confining potential, magnetic 
field, wave vector, frequency, temperature and background dielectric constant. We then varied the period d and 
calculated the determinant det[CV,-(q, ui)] obtained from Eq. p5[) . The zeros of the determinant correspond to the 
collective magnetoplasmon excitations. As is shown in Fig. [2l allowed plasmon resonances appear in the system at 
specific values of the dot spacing d corresponding to the condition det[CV, (q, uj)] = 0. The results of these calculations 
were obtained by taking into account only transitions between the Landau levels n = —1,0,1 . The interdot separations 
where the plasmons can be excited are given by d < 100 A. The excitation spectrum of the collective modes preserves 
the periodicity of the lattice, even with the mixing of the quantum numbers. 



IV. DISCUSSION 



In Fig. [31 we present the dispersion relation for plasmons in a 2D array of quantum dots in a perpendicular magnetic 
field at zero temperature. The parameters used in the calculation are given in the figure. The frequency is plotted 
as a function of q x with q y = 0. As q x increases in the long wavelength limit, the group velocity is initially small but 
then rapidly becomes negative for a range of wave vector. Then, over a narrow range of wave vector, the frequency 
increases before it again decreases as it approaches the Brillouin zone boundary. The finite frequency of excitation in 
the long wavelength limit is due to the allowed transition between the valence and conduction band. In calculating 
the dispersion, we must solve for the zeros of the (2N + l) 4 x (2N + l) 4 -dimensional det[Cy ;(q, u>)] as a function of 
frequency, where i, j = 0, ±1, ±2, • • • , ±N are the energy level labels in the valence and conduction bands. Thus, when 
we include the 0, ±1 levels only, the determinant has dimension eighty-one. We note that the negative dispersion must 
be due to the effective magnetic field since the group velocity is always positive for a graphene sheet in the absence 
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FIG. 2: det[dj (q, w)] as a function of an array spacing d calculated using Eq. (|15[) . We chose the slope of the confinement Ao = 
magnetic field B = IT, £(, = 2.5 and the temperature T=0 K. The values of spacing d resulting in det[CVj(q, a>)] = 
correspond to the existence of magnetoplasmon excitations at chosen frequency and wave vector. 



of magnetic field 0, 0, [l6| ■ The commensurability effect between the electron cyclotron orbits and the period of the 
quantum dot array can be explained by the oscillations in the plasma dispersion shown in Fig. [3l The interaction 
term decreases as the wave vector increases, thereby resulting in the decrease in plasmon excitation energy as q x 
approaches the Brillouin zone boundary. Since each value of inter-dot separation has its own set of eigenstates, the 
dispersion relation could be drastically changed by choosing a different value of d. We note that a zero-gap graphite 
sheet only excites the inter-band electron-hole excitations at zero temperature. At finite temperature, there will be 
intra-band excitations which induce intra-band plasmon modes. 

In summary, the methodological achievement of the approach presented in this paper is the mapping of the single- 
electron eigenfunctions (Eq. @) and eigenenergies in magnetic field in linear confinement in graphene onto the 
eigenfunctions and eigenenergies of a new set of quasiparticles of quasielectrons and quasiholes, without confinement 
in an effective magnetic field depending on the slope of the confinement. The quasiparticle technique is also well suited 
for calculating the different linear responses (e.g., conductivity, linear paramagnetic susceptibility, optical response, 
etc.) which can be found by computing the average (G r (uj, t)G a {lo, — r, )) involving the retarded and advanced Green's 
functions. Also, the calculation of the electronic heat capacity, which, as usual, is found from the the temperature 
derivative of the mean energy, utilizes the quasiparticle energies. Thus, the scope of application of the quasiparticle 
concept is wide and can be employed in both transport, thermodynamic and optical studies. We showed that the 
dispersion equation yields allowed plasmon excitations for specific inter-dot separation at chosen wave vector and 
frequency. Our dispersion curve shows negative group velocity which is generated by the effective magnetic field 
arising from the external magnetic field and confining potential. 
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